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ABSTRACT 

We present optical and near-IR observations of the host galaxy of GRB 020127, for which we measure 
R—Ks = 6.2 mag. This is only the second GRB host to date, which is classified as an ERO. The spectral 
energy distribution is typical of a dusty starburst galaxy, with a redshift, z w 1.9, a luminosity, L « 5 
L*, and an inferred stellar mass of ~ 10^^ — 10^^ M©, two orders of magnitude more massive than 
typical GRB hosts. A comparison to the z ~ 2 mass-metallicity (M-Z) relation indicates that the host 
mctallicity is about 0.5 Z©. This result shows that at least some GRBs occur in massive, metal-enriched 
galaxies, and that the proposed low mctallicity bias of GRB progenitors is not as severe as previously 
claimed. Instead we propose that the blue colors and sub-L* luminosities of most GRB hosts reflect their 
young starburst populations. This, along with the locally increased fraction of starbursts at lower stellar 
mass, may in fact be the underlying reason for the claimed low metallicity bias of z ^ 0.2 GRB hosts. 
Consequently, GRBs and their host galaxies may serve as a reliable tracer of cosmic star formation, 
particularly at z ^ 1 where the M-Z relation is systematically lower, and the star formation rate is 
dominated by sub-L* galaxies. 

Subject headings: gamma-rays: bursts — galaxies:starburst 



1. INTRODUCTION 

Observations of the host galaxies of long-duration 7-ray 
bursts (GRBs) have led to the general consensus that they 
are faint and blue. This has been interpreted as evidence 
for intense star formation activity, as well as low metallic- 
ity and stellar mass (Fynbo et al. 2003; Le Floc'h et al. 
2003; Christensen et al. 2004; Fruchter et al. 2006; Stanek 
et al. 2006). In particular, a comparison of z ^ 0.2 GRB 
hosts to Sloan Digital Sky Survey galaxies suggests that 
GRBs occur preferentially in low metallicity galaxies with 
Z - 0.1-0.5 Z© (Stanek et al. 2006). These authors 
also propose an upper metallicity cutoff of about 0.15 Zq 
for typical GRBs (i.e., with E ~ 10^^ erg). At higher 
redshifts, z ~ 0.5 — 1, a comparison to the host galax- 
ies of GOODS supernovae indicates that GRBs tend to 
occur in lower luminosity galaxies (by about 1 mag) and 
therefore presumably at lower metallicity (Fruchter et al. 
2006). Metallicities measured from afterglow absorption 
spectra (at z > 2) span a wider range, Z ~ 0.05 — 0.5 
Zq (e.g., Berger et al. 2006b; Prochaska 2006). Theoreti- 
cal studies have also argued for "low" metallicity of GRB 
progenitors (MacFadyen & Woosley 1999), although we 
stress that the definition is not well-quantified, and black 
hole remnants can apparently be formed even above solar 
metallicity (Heger et al. 2003). 

These observations raise two crucial and related ques- 
tions: (i) Is there in fact a low metallicity bias for GRB 



progenitors, in both an absolute sense and relative to the 
average metallicity at any redshift? (ii) Can GRBs and 
their host galaxies be used as a representative tracer of star 
formation at high redshift (z ^ 1)? The answer to these 
questions is of great importance given the unique ability of 
GRBs to probe the interstellar medium and star forming 
environments of galaxies over a wide redshift range, ex- 
tending to z > 6 (Berger et al. 2006a). It is also crucial to 
understand the effect of metallicity in light of the redshift 
evolution of the average metallicity and mass function of 
galaxies (e.g., Erb et al. 2006). Equally important, the 
metallicity range for GRB progenitors impacts our under- 
standing of the progenitor population and GRB formation 
scenarios (MacFadyen & Woosley 1999; Heger et al. 2003). 

The prevalence of GRBs in blue galaxies may also reflect 
an observational bias against dusty galaxies due to obscu- 
ration of the optical afterglow. Two lines of evidence sug- 
gest that this may not be a significant problem. First, the 
host galaxies of the "optically dark" GRBs (those with 
clear evidence for dust obscuration in the optical band) 
are typically not redder than the hosts of optically-bright 
GRBs (Berger et al. 2003; Le Floc'h et al. 2003). A notable 
exception is GRB 030115, with a dust-reddened afterglow 
and a host galaxy with R—K k, 5.4 mag (Levan et al. 2006). 
Second, radio, submillimeter, and IR observations have 
not led to a preferential detection of the dark GRB hosts 
(Barnard et al. 2003; Berger et al. 2003; Le Floc'h et al. 
2006), despite an expectation that these galaxies should 
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have obscured star formation. Conversely, the host galax- 
ies with long-wavelength detections (and hence obscured 
star formation), have blue colors typical of the GRB host 
population as a whole (Berger et al. 2003). 

Here we present optical and NIR observations of 
GRB 020127, which reveal that the host galaxy is an ex- 
tremely red object (ERO) with R — Kg = 6.2 mag. The 
GRB position is obtained from X-ray and radio observa- 
tions of the afterglow. The host SED is best fit by an 
obscured starburst galaxy template at z « 1.9, with a 
stellar mass^ of ~ 10^^'^ Mq and an inferred metallicity 
(Erb et al. 2006) of about 0.5 Zg, suggesting that GRB 
progenitors can in fact exist at near-solar metallicity, and 
may occur in the most massive galaxies at z > 1. 

2. AFTERGLOW AND HOST GALAXY OBSERVATIONS 

GRB 020127 was discovered by the High Energy Tran- 
sient Explorer //satellite on 2002 Jan. 27.875 UT, with a 
positional accuracy of 8' radius (Ricker et al. 2002). The 
burst duration and fluence in the 2—400 keV band are 18 
s and 2.7 x lO"*" erg cm~^ (Sakamoto et al. 2005). 

Optical observations did not uncover an afterglow candi- 
date, with R,I> 19.5 mag at 4.4 and 8.2 hr after the burst, 
respectively (Lamb et al. 2002), and R> 21.5 mag at 3.1 
hr over 75% of the error circle (Castro Ceron et al. 2002). 
The Galactic extinction in the direction of GRB 020127 is 
low, E{B -V)= 0.048 mag (Schlegel et al. 1998). 

2.1. X-ray, Radio, and Optical/NIR Imaging 

We initiated two 10 ks observations with the Chan- 
dra X-ray Observatory, 4.14 and 14.64 d after the burst, 
to identify the fading X-ray afterglow. The data were 
obtained with the Advanced CCD Imaging Spectrome- 
ter (ACIS) and reduced and analyzed using the CIAO 
software package^. A comparison of the two epochs re- 
veals three fading sources not associated with bright stel- 
lar counterparts. Of these, only one is detected in both 
epochs, and has faded with high significance, with 2.6 and 
1.0 count ks~^ (0.3 — 7 keV), respectively. Using a Galactic 
column of 3 X 10^° cm~^ (Schlegel et al. 1998) and a pho- 
ton index of 1.5 ±0.4, we find fluxes of (5.5 ±0.1) x 10"" 
and (2.1 ± 0.7) x lO"^** erg cm~^ s""'^, respectively, corre- 
sponding to a power law decay rate, Fi, oc t~°'^, typical of 
GRB afterglows. The position of the X-ray afterglow can- 
didate is a=08*'15'"01.42^ (5=+36°46'33.9" (J2000), with 
an uncertainty of about 1" in each coordinate. 

We observed this source with the Very Large Array 
(VLA9) at a frequency of 8.46 GHz on 2002 Feb. 14.20, 
16.23, 21.97, and Mar. 18.10 UT. The data were processed 
using AIPS. In the first observation we detect an object 
with a flux of 222 ± 63 /xJy, coincident with the X-ray 
position, at a=08'>15™01.42^ (5=+36°46'33.45" (J2000), 
±0.03" in each coordinate. Subsequent observations re- 
veal that the object has faded, with 3a limits of 150 /jxJy 
(Feb. 16.23 and 21.97) and 65 /iJy (Mar. 18.10), suggesting 
that this is the radio afterglow of GRB 020127. 

^We use Vega magnitudes (unless otherwise noted), and the standard cosmology with Ho = 71 km Mpc~^, flm = 0.27, and fljy = 0.73. 
®http: / / cxc. harvard. cdu/ciao 

^The VLA is operated by the National Radio Astronomy Observatory, a facility of the National Science Foundation operated under coop- 
erative agreement by Associated Universities, Inc. 



Optical observations were obtained with the Large For- 
mat Camera (LFC) on the 200-inch telescope at the 
Mt. Palomar Observatory on 2002 Feb. 4.29 and Feb. 6.34 
UT for a total of 4300 s {g'), 2400 s (r'), and 1800 s {i'). 
The images were bias-subtracted, flat-fielded, and coad- 
ded using IRAF, and photometry was performed relative 
to SDSS. At the position of the X-ray and radio after- 
glows we detect a faint object with i' = 23.89 ± 0.13 mag, 
r' > 23.5 mag, and g' > 25.8 mag; limits are 2a and 
magnitudes are in the AB system. We identify this object 
as the host galaxy of GRB 020127. Further observations 
were obtained with the Echellete Spectrograph and Imager 
(ESI) on the Keck II 10-m telescope on 2002 Mar. 13.28 
UT in R (1500 s) and on 2003 Feb. 28.36 UT in I (600 s). 
The data were reduced as described above. We measure 
/ = 23.56 ± 0.10 mag and R = 24.73 ± 0.15 mag. 

NIR observations in Kg and J were obtained with the 
Near Infra-Red Camera (NIRC) on the Keck I 10-m tele- 
scope on 2002 Oct. 13.63 and Nov. 17.48 UT, respectively, 
for a total of 1080 s in each filter. The individual frames 
were dark-subtracted, flat-fielded, and corrected for bad 
pixels and cosmic rays using custom IRAF routines. Pho- 
tometry was performed relative to the standard star Feige 
16. The host galaxy has J = 20.37 ± 0.10 mag and 
Kg = 18.54 ± 0.05 mag. Optical/NIR images of the host 
are shown in Figure 1. 

Finally, on 2002 Apr. 6 UT we observed the host galaxy 
with the Hubble Space Telescope using the Space Telescope 
Imaging Spectrograph (STIS) as part of program GO 9180 
(PI: Kulkarni). A total of 4868 s were obtained with the 
CL filter. We processed and combined the individual ex- 
posures using the IRAF task drizzle (Fruchter & Hook 
2002), with pixfrac=0.8 and pixscale=0 . 5. At the po- 
sition of the afterglow we detect an extended object with 
an AB magnitude of 24.8 ± 0.1 mag (Figure 1). 



2.2. Optical Spectroscopy 

We obtained optical spectra of the host galaxy using 
ESI on four separate occasions. The data were reduced us- 
ing custom IRAF routines to bias-subtract, flat-field, and 
rectify the ten individual echelle orders. Sky subtraction 
was performed using the method and software described 
in Kelson (2003). Wavelength calibration was performed 
using CuAr and HgNeXe arc lamps and air-to- vacuum and 
heliocentric corrections were applied. The spectrum cov- 
ers the range 0.39 to 1.05 /J,m at a resolution of 11.5 km 
s""'^. We detect continuum emission beyond 0.6 /im, but 
no emission or absorption lines are clearly identified. 

A 5400 s spectrum of the host was obtained with the Low 
Resolution Imaging Spectrometer (LRIS) on the Keck 1 10- 
m telescope on 2004 Apr. 22 UT. The wavelength coverage 
is 0.36 to 0.96 /xm with a resolution of 3.3 A (0.36 — 0.58 
^m) and 5.6 A (0.58 — 0.96 /zm). The data were reduced 
as described above. As in the ESI spectra, we detect the 
host continuum but no emission or absorption features. 
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3. HOST GALAXY PROPERTIES 

The optical/NIR spectral energy distribution of the host 
galaxy is shown in Figure 1. With R—Kg = 6.2 ± 0.2 mag 
the host is classified as an ERO, and is the reddest GRB 
host to date. EROs fall in two general categories of old, 
passively-evolving elliptical galaxies, and dust-obscured 
star forming galaxies (McCarthy 2004). Long GRBs have 
never been localized to elliptical galaxies, and we further 
distinguish the two possibilities using the observed SED. 
We use the hyperz package (Bolzonella et al. 2000) to fit 
a range of model templates (starburst, E, SO-Sc, and Irr) 
with the redshift, stellar population age, and extinction 
as free parameters. We assume a Calzetti et al. (2000) ex- 
tinction curve. The elliptical galaxy model provides an un- 
satisfactory fit with a probability of only 0.7% {Xr =4.1; 
Figures 1 and 2), due to the expected flat slope (in F\) 
beyond 1 fxm and a brighter g'-band magnitude. 

The starburst template, on the other hand, provides 
an excellent fit (P^ax = 99.9%; Figure 2). The best-fit 
parameters are z = 1.9^° 4, an age of about 0.7 Gyr, a 
rest-frame extinction. Ay ~ 0.5 mag, and an absolute 
rest-frame magnitude, Mab{B) = —23.5 ±0.1 mag. This 
corresponds to L » 5 L* (Dahlen et al. 2005; Willmer 
et al. 2006). A comparison to UV-selected galaxies at 
z ~ 2 (Shapley et al. 2005) indicates that the host mass is 
M* ~ 10"-i2 Mq. Using the relation of Kcnnicutt (1998) 
and the observed i?-band (2200 A rest-frame) flux of 0.46 
^ Jy, we estimate an unobscured star formation rate of ~ 6 
Mq yr^^. The value corrected for extinction is about an 
order of magnitude larger. 

The host of GRB 020127 is distinguished from other 
GRB host galaxies in several ways. First, it is more than 3 
magnitudes redder than the average value for GRB hosts, 
and even a magnitude redder than the host of GRB 030115 
(Figure 3). This is due to the combined effect of extinc- 
tion and a more evolved stellar population compared to 
the typical value of about 0.1 Gyr for other GRB hosts 
(Christensen et al. 2004). Second, it is four magnitude 
more luminous than the median value of Mab{B) « —19.5 
mag for GRB hosts (with a range of —16 to —22 mag; Fig- 
ure 4). Third, the inferred stellar mass is nearly two orders 
of magnitudes larger than the median value of 10^'^ Mq 
for GRB hosts (Christensen et al. 2004; Savaglio et al. 
2006). Finally, the unobscured specific star formation rate 
of about 1 Mq yr~^ {L/L*)~^, is nearly an order of mag- 
nitude lower than for other GRB hosts (Christensen et al. 
2004), but we note that with the extinction correction it 
is in fact similar. 

Perhaps most importantly, a comparison to the mass- 
metallicity relation of UV-selected galaxies at 2; ~ 2 (Erb 
et al. 2006) indicates that the host of GRB 020127 has 
12 + log (0/H) « 8.6, or about 0.5 Zq. This is also simi- 
lar to the metallicities of the so-called distant red galaxies 
(DRGs; J — K > 2.3 mag), which have median ages and 
masses of ^ 2 Gyr and ^ 10^^ Mq , respectively (Forster 
Schreiber et al. 2004; van Dokkum et al. 2004). 

4. DISCUSSION 

The host galaxy of GRB 020127, identified through X- 
ray and radio afterglow observations, is the reddest, and 



most luminous and massive GRB host discovered to date. 
Unfortunately, the early optical limits (§2) are too shallow 
to ascertain whether the afterglow was significantly extin- 
guished by dust; from the X-ray flux at 4.14 d we expect 
i? ~ 21 mag at the time of the early optical observations, 
at the level of the available limits. For GRB 030115, the 
only other burst with an ERO host galaxy, the afterglow 
itself was shown to be significantly dust extinguished, with 
R—K « 6 mag (Levan et al. 2006). It is therefore possible 
that dusty galaxies are under-represented in the current 
GRB host sample, due to obscuration of the optical after- 
glow, but we stress that the hosts of several dark bursts 
are not significantly redder than the median of the sample 
(Berger et al. 2003). 

The discovery of a GRB in a massive, and hence metal- 
enriched, starburst galaxy indicates that the proposed low 
metallicity bias of GRB progenitors (^ 0.1 Zq; Pruchter 
ct al. 2006; Stanek ct al. 2006) is not likely to be as severe 
as previously claimed, particularly at z ^ 1. This appears 
to be true not only in terms of an absolute metallicity 
cutoff, but also relative to the mass-metallicity relation of 
galaxies at high redshift. In particular, while the hosts 
of GRBs at z ^ 0.2 appear to have metallicities at the 
low end of the distribution for local galaxies (Stanek et al. 
2006), at least some GRBs at 2; ~ 2 occur in the most 
metal-enriched galaxies at that redshift. This conclusion 
is also supported by the metallicities derived from after- 
glow absorption spectra (at z > 2), which range up to 
solar values (Berger et al. 2006b; Prochaska 2006), as well 
as by the detection of submillimeter emission from some 
GRB hosts (Berger et al. 2003). Since GRB progenitors 
appear to occur at least up to ~ 0.5 Zq, and given that 
the M-Z relations at z ~ 1 and 2 are systematically 
lower by about 0.3 and 0.6 dex, respectively, compared 
to the local relation (Tremonti et al. 2004; Savaglio et al. 
2005; Erb et al. 2006), we conclude that even if a slight low 
metallicity bias does exist, its effect will diminish beyond 
z ~ 1. 

Instead, the blue colors of GRB host galaxies may re- 
flect their young stellar populations. This is likely related 
to the fact that GRB progenitors are massive stars, which 
explode within a few million years of formation, thereby 
leading to the selection of young starburst galaxies. In 
fact, for z ^ 1 — 2, the average R — K color of a 0.1 Gyr 
population is about 2.5 mag, in good agreement with the 
observed colors of GRB hosts, while for 0.3 and 1 Gyr 
populations it is about 1 and 2.5 magnitudes redder, re- 
spectively (Figure 3). 

This explanation might also account for the transition at 
low redshift (z <^ 0.2) from a representative to a predomi- 
nantly low luminosity host population, since locally young 
starburst activity occurs primarily in low mass galaxies. 
Brinchmann et al. (2004) show that locally the fraction of 
galaxies with recent starburst activity increases strongly 
with decreased mass. Similarly, Bell et al. (2005) show 
that while at z ~ 0.7 nearly half of all galaxies with 
M ^ lO"'^" Mq undergo intense star formation activity, 
locally, this number is less than 1%. In this scenario, the 
driving parameter is stellar population age, which may be 
misidentified as a low metallicity bias. Clearly, a study of 
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the stellar population ages of the low redshift GRB hosts 
is required to assess whether age, and not metallicity, in 
fact plays the underlying role. 

Thus, we conclude that at least at ^ 1, GRBs and 
their host galaxies are likely to trace star formation in 
a relatively unbiased way, with possibly a preference for 
younger starburst populations. At z ^ there seems to 
be a bias in favor of low luminosity /mass galaxies (Stanek 
et al. 2006), but this may be driven by stellar population 
age and not metallicity. The luminosity function of GRB 
hosts (see also Jakobsson et al. 2005) therefore indicates 
that the bulk of the star formation at 2: ^ 1 takes place in 
galaxies fainter than L*, in good agreement with the steep 
faint-end slope of the luminosity function of high redshift 



galaxies. 
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Fig. 1. — Spectral energy distribution of the host galaxy of GRB 020127 along with the best fit starburst and elliptical 
galaxy models. The combination of a blue J — color and a deep upper limit in the 5-band indicates that the starburst 
template is a much better fit compared to an elliptical galaxy. The top panel shows the host galaxy in each of the five 
observed bands. 
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Fig. 2. — Probability distribution (black lino) and reduced (gray line) for the host galaxy photometric redshift 
calculated with the hyperz software package (Bolzonella ct al. 2000) using a starburst template. The best fit redshift is 
z = 1.9t.oA- The inset shows the reduced ^or the best fit solution using the various galaxy templates. Clearly, the 
starburst template provides the only adequate fit to the data. 
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Fig. 3. — Apparent R — Kg color plotted versus rcdshift for GRB host galaxies (filled circles; Chary et al. 2002; Le 
Floc'h et al. 2003; Berger ct al. 2003; Levan et al. 2006), LBGs (open circles; Shapley et al. 2001; Steidel et al. 2004), 
submillimeter galaxies (stars; Small et al. 2004), and GOODS galaxies (shaded; Somerville et al. 2004). Arrows designate 
GRB hosts without a measured redshift. The host galaxy of GRB 020127 (black square) is significantly redder than the 
typical color of R—Kg ~ 2.5 — 3 mag. The gray lines are tracks of R—Kg color as a function of rcdshift from the population 
synthesis code of Bruzual & Chariot (2003) for a stellar population age of 0.1 Gyr (solid), 0.3 Gyr (dashed), and 1 Gyr 
(dot-dashed). The majority of the GRB hosts appear to have ages of 0.1 — 0.3 Gyr. These young stellar population 
ages, and not metallicity, may be the underlying reason for the blue colors and low luminosities of GRB hosts. 
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Fig. 4. Distribution of absolute rcst-franic _B-band magnitudes for GRB host galaxies (Berger et al. in prep.). The 
median value for the sample is Mab(-B) = —19.5 mag, or 0.1 L*, and it does not appear to significantly vary between 
z <1 and z > 1 (values in the figure are the median and standard deviation). GRB 020127 clearly stands out as the most 
luminous, and thus most massive host galaxy detected to date, but we note that the overall luminosity distribution is in 
good agreement with that of other high redshift galaxy sample (see also Jakobsson et al. 2005). The dashed curve is a 
representative Schechter luminosity function with Mg = —21 mag and a = —1.4. 



